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Abstract. Using quasiconformal mappings, we prove that any Rie- 
mann surface of finite connectivity and finite genus is conformally equiv- 
alent to an intrinsic circle domain f2 in a compact Riemann surface S. 
This means that each connected component B of S \ Q is either a point 
or a closed geometric disc with respect to the complete constant curva- 
ture conformal metric of the Riemann surface (J! U B). Moreover the 
pair (fi, S) is unique up to conformal isomorphisms. We give a gener- 
alization to countably infinite connectivity. Finally we show how one 
can compute numerical approximations to intrinsic circle domains using 
circle packings and conformal welding. 



1. Introduction 

Let Q be a finitely connected domain in the Riemann sphere C. A classical 
theorem of Koebe states that O is conformally equivalent to the complement 
of a finite set of pairwise disjoint closed discs and points. Such a domain is 
called a finitely connected circle domain. Koebe's theorem was extended by 
He and Schramm [3] to apply to domains with countably many complemen- 
tary components. Schramm later gave a different proof of this result using 
transboundary extremal length [TJ. 

Theorem A. [H Theorem 0.1] Let Q be a domain in C such that boundary 
d£l has at most countably many components. Then 0, is conformally home- 
omorphic to a circle domain Q* in C. Moreover, 0* is unique up to Mobius 
transformations and every conformal automorphism of is the restriction 
of a Mobius transformation. 

If SI is a domain in C that has at least three complementary components, 
or at least one complementary component that is not a puncture, then we 
say that f2 is hyperbolic. The reason is that there is then an unbranched 
analytic covering map from the unit disc D onto £1, which can be used to 
transfer the Poincare metric on D to a complete conformal metric of constant 
curvature —1 on 0, which is called the hyperbolic metric. 

Every ring domain is conformally equivalent to a round annulus, a punc- 
tured disc, or the punctured plane. For a round annulus C \ (B\ U B2), each 
complementary component B{ is a spherical disc, i.e. a closed ball in the 
spherical metric. Observe that B{ is also a closed ball with respect to the 
hyperbolic metric belonging to the domain QuBi, which is a larger spherical 
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disc. In this paper we generalize this property to obtain a new canonical 
form for multiply connected domains. 

Theorem 1. Let Vt be a finitely connected domain in the Riemann sphere 
C with complementary components K\, . . . , K n . Suppose that each domain 
Q U Ki is hyperbolic. Then VL is conformally equivalent to a domain £1* in C 
with complementary components L±, . . . ,L n such that for each i = 1, . . . , n, 
0, U Li is hyperbolic and either Li is a puncture or Li is a closed disc with 
respect to the hyperbolic metric o/OULj. Moreover, f2* is unique up to 
Mobius transformations. 

We call the canonical domains given by this theorem intrinsic circle domains. 

Theorem [1] has a generalization to positive genus ambient surfaces in 
place of C, which we will now explain. In any Riemann surface S, a closed 
geometric disc will mean a closed ball of some positive radius with respect 
to the appropriate complete conformal metric of constant curvature on S, 
(i.e. hyperbolic, Euclidean or spherical), with the extra condition that it 
must be homeomorphic to the closed unit disc. Equivalently, the radius of 
the ball must be strictly less than the injectivity radius of the metric at the 
center of the disc. A circle domain Q* in S is a connected open subset of 
S for which each complementary component is either a point or a closed 
geometric disc in S. Via the uniformization theorem, He and Schramm 
extended Theorem [A] to this setting, as follows. 

Theorem B. [U Theorem 0.2] Let O be an open Riemann surface with finite 
genus and at most countably many ends. Then there is a closed Riemann 
surface R* such that Q is conformally homeomorphic to a circle domain Q* 
in R* . Moreover the pair (R*,£l*) is unique up to conformal homeomor- 
phisms. 

Likewise, we can extend Theorem[T]to deal with the case of arbitrary finite 
genus. For the moment we also restrict ourselves to finite connectivity; we 
will relax this condition in £j5j 

Theorem 2. Let Q be a Riemann surface of finite genus and finite connec- 
tivity. Then there is a conformal embedding (poffl into a compact Riemann 
surface S, of the same genus as 0, so that S \ </?(f2) is the union of dis- 
joint closed, connected and simply connected sets L\, . . . , L n , and for each 
i 6 {1, . . . , n}, Li is either a single point or a closed geometric disc with re- 
spect to the Riemann surface tp(£l) ULi. Moreover, the pair (S,ip) is unique 
up to conformal homeomorphisms. 

The condition that Q and S have the same genus means that S has the 
minimal possible genus among all compact Riemann surfaces into which Q 
may be embedded. We will refer to the domain <p(£l) as an intrinsic circle 
domain in S. Note that Theorem [2] includes Theorem Q] as the special case 
of genus 0, so this terminology is consistent; we shall only give a proof of 
Theorem [2j 

There are some simple special cases. For example, there is precisely one 
case in which the appropriate geometry of (p(£l) U Lj is spherical; this oc- 
curs when $7 is simply connected, in which case the statement reduces to 
the Riemann mapping theorem. The doubly connected genus cases occur 
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when f2 is conformally equivalent to the punctured plane C*, the punctured 
disc B*, or a round annulus. Note that in the case of the punctured disc, 
the appropriate geometry for (p(£l) U Li is Euclidean for one complemen- 
tary component and hyperbolic for the other. Another case in which more 
than one type of geometry must be considered occurs when Q is a simply 
connected domain with two points removed; then the resulting canonical 
domain is either a triply-punctured sphere or it is C* \ K, where K is the 
image under the exponential map of a closed disc of radius less than ir. If 
f2 is a Riemann surface of genus 1 with one end, then f2* is the comple- 
ment of a point or of a Euclidean disc in C/A, where A is a lattice in C. 
In all other cases the natural geometry of (p(£l) U L.; is hyperbolic for every 
complementary component Lj. 

Our definition of intrinsic circle domains was motivated by an observa- 
tion about extremal cases in the Grotzsch or Polya-Chebotarev problem, 
which asks for the minimizer of the logarithmic capacity among compact 
connected sets containing a given finite set of points in C. This problem 
often appears in the course of studying other extremal problems in geo- 
metric function theory. For example, there are recent applications to the 
Bloch-Landau constant [T] and to Smale's mean value conjecture [2]. It is a 
classical result of Lavrentiev and Goluzin that the extremal continuum E is 
unique and is the union of finitely many analytic arcs A4, which are trajec- 
tories of a certain rational quadratic differential (see for example [3])- It is 
not hard to show that each arc Ai is a geodesic arc in the hyperbolic metric 
of the domain (C\E)L)Ai. A similar condition arises for local minima of the 
condenser capacity for domains separating one finite set of points from an- 
other. In [5] it is observed that the extremal continuum E enjoys harmonic 
symmetry. This means that for any subarc I of E, the harmonic measure of 
the complementary domain C\E with respect to the point 00 assigns equal 
masses to each side of /. Note that both the geodesic arc condition and the 
harmonic symmetry condition have the property that they may be verified 
by checking each arc Ai separately. To check Ai, we need to know only the 
conformal class of the pair (C\E, Ai). Our notion of intrinsic circle domains 
arose by analogy with this property. 

We now outline the rest of the paper. Section [2] gives a simple qualitative 
distortion bound for quasiconformal extensions of conformal maps between 
ring domains, which we later use several times. Section [3] gives the proof 
of uniqueness in Theorem [2j and section 2] gives the existence proof. In 
section [5l we extend Theorem [2] to include some cases with countably infin- 
itely many complementary components, subject to a geometric constraint. 
Section [6] deals with a mixed condition, in which some boundary compo- 
nents are required to be circles in the spherical metric while the others are 
required to be intrinsic circles in the sense of Theorem [TJ In section [7] we 
discuss the use of circle packings for the numerical approximation of intrinsic 
circle domains, and illustrate with some examples. 
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2. A GEOMETRIC LEMMA 

Lemma 3. Suppose A\ and A2 are ring domains in C whose inner bound- 
aries are circles C\ and C2 respectively. Suppose there is a conformal home- 
omorphism F : A\ — > A2 under which C\ corresponds to C2. Then the 
induced homeomorphism F : C\ — > C2 has a K -quasiconformal extension 
between the interior discs of C\ and C2, where the constant K depends only 
on the conformal modulus of A\ . 

Proof. Lemma [3] is a consequence of [U Thm 1.4], but it can also be proved 
simply as follows. After applying similarities, we may assume that C\ and C2 
are both the unit circle <9D. Let m be the conformal modulus of A\ . Apply 
Schwarz reflection across <9B to extend F to a conformal homeomorphism 
F between ring domains A\ and A%. If m = 00 then F is a conformal 
homeomorphism from the punctured plane to itself, so is in fact a similarity, 
and its restriction to 9B is an isometry. Otherwise, the doubled ring domains 
A{ are hyperbolic, F is a hyperbolic isometry, and the hyperbolic length of 
c© in each domain depends only on m. We claim that the density of the 
hyperbolic metric of A\ on <9B is bounded above and below in terms of m. 
Indeed, the hyperbolic length of <9B is a function of m so we get such a 
bound by applying the Koebe distortion theorem to a single-valued lift of 
logF to the universal cover of A\. 

It follows that F : <9IB — > dO is bi-Lipschitz , with constants that depend 
only on m. Therefore the radial interpolation of the boundary correspon- 
dence provides a quasiconformal extension with dilatation bounded in terms 
of the modulus of A\. □ 

3. Uniqueness 

Suppose that O and 0' are two intrinsic circle domains in closed Riemann 
surfaces S and S' respectively. Let / : f2 — > £1' be a conformal homeomor- 
phism. We have to show that / is the restriction of a conformal isomorphism 
/ : S -> S'. 

First, suppose that O is simply connected or doubly connected. If £1 is 
simply connected or is doubly connected with infinite conformal modulus, 
then S is the Riemann sphere and £1 is a circle domain. 0' has the same 
connectivity as f2 since they are homeomorphic, so 0,' is also a circle domain 
in the Riemann sphere S. The uniqueness part of Koebe's theorem then 
shows that / is the restriction of a Mobius map. Otherwise, f2 is a ring 
domain of finite modulus, as is O'. Then S and S' are both of genus and 
we have to show that / extends to a Mobius map. Consider a complemen- 
tary component B\ of f2, and let B[ be the corresponding complementary 
component of O'. Let ir : D — > 0, U B\ and it 1 : D — > £1' U B[ be any Riemann 
maps. Then by hypothesis 7r _1 (i?i) and tt'~ 1 (B' 1 ) are closed hyperbolic discs 
in D, which are bounded by Euclidean circles. Thus 7r _1 (f2) and 7r' _1 (Q / ) are 
round annuli. Any conformal homeomorphism between round annuli is the 
restriction of a Mobius map, so we can extend 7r /_1 o / o tt to a Mobius map 
g. Now define f\ = n' o g o 7r _1 . Then f\ : U B\ — > £2' U B[ is a conformal 
homeomorphism extending /. Similarly, / extends to a conformal homeo- 
morphism f2 : £1 U B2 — > B' 2 . Gluing f± and /2 together we obtain a 
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conformal homeomorphism from the Riemann sphere to itself which extends 
/. This must be a Mobius map. Since is conformally equivalent to some 
round annulus, which is an intrinsic circle domain, this in fact shows that 
is a round annulus. 

Now suppose that is at least triply connected. Consider any comple- 
mentary component Bi of fi, and let B[ be the corresponding complementary 
component of 0'. This makes sense since the homeomorphism / induces a 
bijection between the ends of Vt and the ends of 0,'; for an intrinsic circle 
domain each end corresponds to precisely one complementary component 
since the complementary components are all contractible. 

Since f2 and f2' are at least triply connected, U and fi' U B[ are 
not simply connected. Therefore we must consider their universal covers, 
U and W respectively, in order to understand their hyperbolic metrics. Let 
7r : U —7- QUBi and it' : W — > Q'UB^ be unbranched analytic covering maps, 
with deck transformation groups Tj and T\ respectively. Since Bi and B[ are 
contractible subsets of f2 U Bi and Q' U B[ respectively, the homeomorphism 
/ induces a homotopy equivalence f2 U Bi — > Q' U B^, which in turn induces 
an isomorphism pi : Tj — > T^. Then / lifts to a conformal homeomorphism 

f l :U\7r- 1 (B l )^U'\7r'- 1 (BQ . 

Note that fi is (Tj, r£)-equivariant: for any element 7 £ Tj we have 

Piin) ° fi = fi°l ■ 

The connected components of 7r _1 (Bi) and of ir'^ 1 (B^) are disjoint closed 
discs because and f2' are intrinsic circle domains. Therefore fi is a con- 
formal homeomorphism between circle domains. Hence by the uniqueness 
part of Theorem [B] it is the restriction of a Mobius map Mj . We find that 
Mi takes U onto W since each of the two circle domains has only one non- 
isolated boundary component. In particular, U = W . Moreover, the map 
Mi :U —>W is also (r^, r£)-equivariant, since Pi(^y) o Mi o 7 -1 : U — y U' is a 
Mobius map extending fi and is therefore equal to Mj. It follows that Mi 
descends to a conformal homeomorphism /, : U Bi — > 0,' U B[ that extends 
/■ 

Gluing the extensions fi together for i = 1, . . . , n, we obtain the desired 
conformal homeomorphism / : S — > S' extending /. This completes the 
proof of uniqueness. 

4. Existence 

The first step is to apply Theorem [Bi to map via a conformal homeo- 
morphism onto a circle domain $7* in some compact Riemann surface R. 

Our goal is to construct a quasiconformal homeomorphism / of R onto 
another Riemann surface S so that / is conformal on 0* and the image 
/(O*) is an intrinsic circle domain in S. In fact we will construct a Beltrami 
coefficient /i on R with H^Hoo < 1; then by the measurable Riemann mapping 
theorem we will obtain a Riemann surface S and a quasiconformal homeo- 
morphism / : R — )• 5 such that p(z)f z = fz a.e. on R. By construction, fi 
will be identically zero on Q* , so that / is conformal there. 
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Figure 1. Domains, maps and group actions in the exis- 
tence proof, illustrating the case where R is hyperbolic. 

Consider any connected component B of R\£l* that is not a single point. 
Our aim is to construct \x on R so that f(B) will be a closed geometric disc 
with respect to the Riemann surface U B). Since this is a condition on 
the conformal structure of /(fi* U B), it depends only on the restriction of 
[i to J)* U B, which by construction will be non-zero only on B. This is a 
key point, for it means that we can correct the conformal structure on each 
complementary component separately and the corrections will not interfere 
with each other. This is in contrast to the proof of Koebe's theorem by 
iterated Riemann mapping, where each complementary component has to be 
corrected infinitely many times and the required mapping is obtained in the 
limit. By hypothesis Q* has only finitely many complementary components 
in R, so to ensure that ||^||oo < 1 it will suffice that \\h\\b < 1 for each 
complementary component B. 

Let 7r : U — > 0* U B be an unbranched analytic covering map, where U 
is one of D, C, or C. We rule out the case U = C, since in this case Q* 
is an open disc in C, so it is already an intrinsic circle domain. Let T be 
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the deck transformation group of tt, i.e. the (infinite) group of conformal 
automorphisms 7 : U — > U such that tt o 7 = n. The preimage tt~ 1 (B) has 
infinitely many connected components, each of which is a topological closed 
disc, bounded by an analytic Jordan curve. 

Apply Theorem lAl to the domain 7r _1 (f2*) = U\ir~ 1 (B). This provides a 
conformal mapping h : 7r~ 1 (J7*) — > V C C, where V is a circle domain. Each 
component of 7r _1 (i?) corresponds under h to a component of C\ V that is 
a closed disc, not a singleton, since it is isolated and cannot be separated 
from the other components by a ring domain of arbitrarily large modulus 
contained in tt" 1 (0,*). The remaining end of 7r _1 (Q*) may correspond to 
either a point or a disc in the complement of V; this component is distin- 
guished because it is not isolated. Let V be the union of V and all of its 
isolated complementary components. 

For any non-trivial element 7 € T, the map h o 7 is also a conformal 
mapping of 7T" 1 (fT) onto the circle domain V, so by the uniqueness part of 
Theorem [A] we have /107 = 7'o/?, for a unique Mobius map 7', which restricts 
to an automorphism of V with no fixed points since 7 has no fixed points 
in U. The map that sends 7 to 7' is therefore an injective homomorphism 
T — > T', where T' is a subgroup of Aut(V). Since F' acts simply transitively 
on the components of /i(7r _1 (i?)) and has no fixed points in V, it acts freely 
on V. Since Aut(V) acts properly discontinuously on V, and each compact 
subset of V intersects only finitely many components of h(ir~ 1 (B)), we also 
conclude that T' acts properly discontinuously on V. 

Let C be a connected component oi U \ 7r _1 (O*), and let C" be the 
corresponding component of V \ V, so that C is bounded by a circle dC. 
Let 7Tr : W — > R be an analytic universal covering of R, so that 7r^ (f2*) is 
a circle domain. Let B be a connected component of iT^ i 1 (B), bounded by 
a circle dB. 

We will examine the behavior of the maps h, tt and 7Tr in the neighborhood 
of the curves dC and dB. For the present proof the aim is merely to show 
that h may be extended continuously to a quasiconformal homeomorphism 
h! :U —¥ V that is equivariant with respect to T and V. In fact we will work 
slightly harder, using Lemma [3] to show that the quasiconformal dilatation 
of this extension can be made to depend only on a conformal invariant of 
0. This will be useful in the next section. 

Let a be the simple closed geodesic in the hyperbolic metric of 0* that 
separates B from all the other components of R\Q*. Let E be the ring 
domain bounded between a and B. It has modulus mod(E) = 7r 2 /2£(a), 
where £(a) is the length of a in the hyperbolic metric of 0,*; note that this 
length is a conformal invariant of $7. 

The restriction of ttr to this circle domain is an unbranched analytic 
covering of 0* and is therefore an isometry from the hyperbolic metric of 
w^ 1 (il*) to the hyperbolic metric of Thus the connected components of 
7r^ 1 (a) are simple closed geodesies in the hyperbolic metric of ttJ^ 1 (SI*) and 
are therefore disjoint, each separating a connected component of 7r^ 1 (S) 
from all of the other components of W \ vr^ 1 (^*). It follows that ttr maps 
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each connected component of ir R (E) bijectively onto E. One of these is a 
ring domain E surrounding B. 

Likewise, the restriction of iroh^ 1 to V is an unbranched analytic covering 
of f2*, so 7r o h^ 1 maps each connected component of h o ir~ 1 (E) bijectively 
onto E. One of these connected components is a ring domain E' surrounding 

a. 

Consider the branch of 71"^ 1 o tt o h- 1 that maps E' to E. It is a home- 
omorphism of ring domains taking the inner boundary circle dC onto the 
inner boundary circle dB. Lemma [3] gives an extension to a quasiconformal 
homeomorphism Q from E' U C to E U B. 

Now Q _1 o7r^ 1 gives us a quasiconformal homeomorphism g : B — > C that 
continuously extends the boundary correspondence induced by of h o ir^ 1 . 
We define h\b to be the Beltrami coefficient gz/gz- 

When we solve the Beltrami equation to obtain / : R — > S, the covering 
map n' : V ->• (S\ f (fi*)) U /(B) such that / o vr = vr' o /» will be an 
unbranched analytic covering, mapping the disc C onto f(B), as required. 

5. The case of countably infinite connectivity 

It is natural to ask whether the notion of intrinsic circle domains can be 
extended to domains of countably infinite connectivity. 

5.1. All ends isolated. First, we deal with the case of a domain Q such 
that every end of R is isolated, i.e. has a neighborhood that meets no other 
end. If is embedded in a compact ambient Riemann surface then this 
implies that has only finitely many ends and finite genus, and we already 
understand this case. 

For a more interesting example, let R be a Z-cover of a compact genus 2 
Riemann surface, so that R has infinite genus and two ends. Then let Q be 
a domain obtained by removing countably many disjoint closed topological 
discs with no accumulation point from R. We could hope to modify the 
structure of R on a neighborhood of each disc in order to make the com- 
plementary components be intrinsic discs. However, there is no way to fix 
the two ends of R so that they are represented by intrinsic discs, for the 
resulting ambient Riemann surface would be compact and therefore have 
only finite genus. 

Suppose that is infinitely connected or has infinite genus. We call an 
end of f2 fixable if it has a neighborhood that is a ring domain. 

Lemma 4. For any Riemann surface £1, there is a conformal embedding of 
f2 into another Riemann surface R such that every connected component of 
R\Q is an intrinsic disc and no end of R is fixable. 

Proof. We can find disjoint open neighborhoods of all the fixable ends, for 
example by cutting along closed geodesies or horocycles around each end. 

In the existence proof for Theorem [2] we only used a local surgery to mod- 
ify the ambient Riemann surface in a neighborhood of each complementary 
component. The initial step of passing to a circle domain was technically 
convenient (and will be needed later), but was not really necessary. All we 
really needed to know was that for each individual end E, we can embed f2 
in a Riemann surface Re such that Re\Q is a geometric disc in Re- This is 
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true for any fixable end E, since it has a neighborhood that is a ring domain 
and therefore it has a neighborhood Ue conformally equivalent to a round 
annulus A = {z € C | a < \z\ < 1}, for some < a < 1. Then we can glue 
the open unit disc to $7, identifying A with Ue, to obtain a Riemann surface 
R' E in which is embedded so that the end E corresponds to a connected 
component K of R' E \ Vt. Then we can apply Theorem VK\ to the preimage 
of in the universal cover of R' E , to find out how to modify the conformal 
structure on K to obtain an embedding of Q in a new Riemann surface Re 
so that the end E corresponds to a connected component of Re \ Q that is 
an intrinsic disc. 

To construct R we glue together all the Riemann surfaces Re correspond- 
ing to fixable ends E by identifying the embedded copies of £1. The resulting 
R is locally a Riemann surface, is connected, and is second countable because 
Q can only have countably many fixable ends. □ 

5.2. Non-isolated ends. We now return to the case of a subdomain SI 
of a compact Riemann surface, but assume that f2 is countably infinitely 
connected. It must now be the case that some complementary components 
are not isolated. There is a potential topological obstruction associated 
with complementary components that are not isolated: they cannot be rep- 
resented by intrinsic discs of finite radius. We might attempt to salvage 
something by allowing complementary components to be horodiscs with re- 
spect to the hyperbolic metric, but this would not help in the case of a 
circle domain 17 C C in which some circular complementary component B 
has precisely two points on its boundary that are accumulation points of 
other complementary circles. 

Therefore in the countably infinitely connected case, we define an intrinsic 
circle domain to be a subdomain of a compact Riemann surface R such 
that any complementary component that is not isolated is a singleton and 
every non-singleton complementary component L is a closed geometric disc 
with respect to f2 U L. This definition is intended to be analogous to the 
definition of a circle domain in the countably connected case. 

In order to obtain a positive theorem, we can place a simple conformal 
geometric constraint on the domain. We will say that a countably connected 
domain in a compact Riemann surface R is uniformly separated if there 
exists e > such that each connected component of R \ Q either is a single 
point or is separated from all the remaining complementary components by 
a ring domain of modulus at least e embedded in R \ Q. 

Lemma 5. Suppose VL is a countably- connected uniformly separated domain 
in a compact Riemann surface R. Let (p be any conformal embedding of S7 
into a compact Riemann surface R' of the same genus as R. Then the non- 
singleton complementary components of R' correspond to the non-singleton 
complementary components of R, so (p(£l) is uniformly separated in R' . 

To explain the significance of this, we first note that a complementary 
component B is isolated if and only if Q contains a ring domain with B 
as one of its complementary components; this is a topological condition 
on 0, so the corresponding complementary component in (p(fl) is also not 
isolated. If is uniformly separated, then the connected components of 
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R\Q that are not singletons must be isolated. However, the remaining 
complementary components need not be punctures; some or all of them could 
be accumulation points of sequences of other complementary components. 
For a general open Riemann surface O of countable connectivity and finite 
genus, it is possible for an end of to be represented in one conformal 
embedding by a complementary component that is a non-isolated singleton, 
yet in some other conformal embedding to be represented by a non-singleton. 




Figure 2. Two conformally equivalent countably connected domains. 

We can construct an example of this cavitation behavioiQ as follows. Con- 
sider a subdomain D of C\B obtained by removing countably many complex 
conjugate pairs of circular arcs, where one arc of the n th pair is 

A n = { (1 + 2~ n ) e u : 6 n < t < ir - 6 n } , 9 n \0. 

The domain 

D n = C\ (bu |J (AjUAj^j 

is conformally equivalent via a conformal map that fixes z = 2 to a slit do- 
main S n , which may be obtained from C\{0} by removing finitely many pairs 
of intervals I n , I n contained in the imaginary axis, arranged symmetrically 
about 0, together with an interval (— 6 n i, 6 n i). See figure [2] for a schematic 
illustration. One Ccin show th.a,t if the angles 6 n are chosen to decrease to 
sufficiently fast then the domains S n converge in the Caratheodory topology 
to a limit domain S in which is a non-isolated singleton complementary 
component, but every other complementary component is an isolated inter- 
val of the imaginary axis. Then D is conformally equivalent to S and we 



The term cavitation refers in fluid dynamics to the sudden formation of a bubble as 
a dissolved gas comes out of solution around a nucleation site. 
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have an example of cavitation: the singleton complementary component {0} 
of S corresponds to the complementary component ID of D. 

The idea is that by making n small we can ensure that the extremal 
length of the family of curves in D joining A n to A n is as small as we like. 
This means we can make the gap between the intervals I n and I n as short 
as we like in comparison to the length of /„. In fact to ensure that the end 
represented in D by the unit circle is represented in S by a singleton, it 
suffices to take 6 n = 2~ n . 

We will now prove Lemma 03 which says that uniform separation prevents 
cavitation. 

Proof. Suppose each non-singleton connected component Bi of R \ Vl is sep- 
arated from the remaining complementary components by a ring domain 
Ai C Q, of modulus e. Then the core curve 7, of the ring domain Ai has 
length 27r 2 /e with respect to the hyperbolic metric of Ai. By the Schwarz- 
Pick lemma, the length of 7, with respect to the hyperbolic metric of fl is 
also at most 2ir 2 /e. Consider the simple closed geodesic gi (with respect to 
the hyperbolic metric of Q) that lies in the free homotopy class of 7$. Since 
gi is a length minimizer in its free homotopy class, it is no longer than 7, 
in the hyperbolic metric of It follows that the ring domain Ei bounded 
between gi and Bi has modulus at least e/4. Our reason for passing from 
the ring domains Ai to the ring domains Ei is that the closures of the Ei 
in R are pairwise disjoint, which was not necessarily true of the Ai, and 
moreover the Ei are intrinsically-defined subdomains of 0. 

We now cut fi along each gi to leave a domain U = Cl\{jEi. We claim 
that for any conformal embedding <p of Q into a compact Riemann surface 
R' of the same genus as R, there exists a quasiconformal homeomorphism 
Q : R — >• R! such that f\u = Q\u- To prove this, consider a component 
Di = Bi U Ei of R \ U. It corresponds to a component D\ = B[ U <p>(Ei) of 
R' \<p(U), where B[ is a component of R! \ ip(Q). We claim that there is a 
-RT-quasiconformal homeomorphism Di — > D[ whose restriction to the curve 
gi agrees with ip. Here K depends only on e. To see this, let tp : Di — > B and 
X '■ D'i — » B be Riemann mappings, and note that dDi and dD[ are analytic 
Jordan curves. Apply Lemma [3] to the conformal homeomorphism 

: Y> {Ei) -> X (El) 

and apply a Schwarz reflection to obtain a i^-quasiconformal homeomor- 
phism h : B — )■ B whose boundary correspondence agrees with that of 
Xo ip o where K depends only on e. Now 

X' 1 o h o V : A -> D| 

is a i^-quasiconformal homeomorphism agreeing with tp on the boundary, 
as required. 

Gluing together ip and the X-quasiconformal extensions to the compo- 
nents Di, we obtain a X-quasiconformal homeomorphism 

H-.nu\jBi -> y>(n)u(jB-. 

Since is countably connected, the complement X = R \ (CI U {J Bi) is 
a countable union of singletons, and therefore is a removable set for the 
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solution of the Beltrami equation. Thus H extends to a homeomorphism 
R — > R' , and hence the complementary components of R' that are not among 
the Bi are singletons. 

□ 

Theorem 6. Suppose is a domain in a compact Riemann surface R such 
that has countably infinite connectivity, has the same genus as R, and 
is uniformly separated. Then there is a conformal embedding (p of in a 
compact Riemann surface S so that </?(0) is an intrinsic circle domain in S. 
The pair (S, <p) is unique up to conformal isomorphism. 

Proof. The existence proof is similar to that of Theorem [2J Lemma [5] shows 
that when we apply Theorem [Bj the resulting countably connected circle 
domain CI* C R satisfies the hypotheses of Theorem [6l in addition to having 
the property that its non-singleton complementary components are geomet- 
ric discs. Then we note that by Lemma [3] the quasiconformal distortion of 
the extension C — > B can be bounded in terms of e. Then the resulting 
Beltrami coefficient [i has ||/u||oo < 1, so we may still apply the measurable 
Riemann mapping theorem to obtain the desired conformal structure on S. 
Then Lemma [S] ensures that the non-isolated complementary components of 
f(Cl*) are singletons, so that f(£l*) is an intrinsic circle domain. 

The uniqueness proof goes through as before, but with one additional 
step. Because of the separation hypothesis there is no difficulty in extending 
a conformal homeomorphism / : Q — > Q' between two countably-connected 
intrinsic circle domains across each non-singleton complementary compo- 
nent. We thus obtain a conformal homeomorphism between domains whose 
complements consist of countably many points, which in turn extends to 
a conformal homeomorphism / : S — >■ S' between the ambient Riemann 
surfaces. □ 

It is likely that the condition of uniform separation can be weakened 
by making use of results on existence and uniqueness of solutions to the 
Beltrami equation in the case where the Beltrami coefficient has norm 1 but 
the measure of the set on which the distortion is large is controlled. 

6. Mixtures of intrinsic and extrinsic circles 

Now we present a generalization of Theorems [2] and [6] in which the con- 
ditions on the complementary components of O in S are mixed. That is, 
some components are punctures or geometric discs in the natural geometry 
of the ambient compact Riemann surface S, while each of the others is an 
intrinsic disc, i.e. a geometric disc with respect to its own union with f2. 

Theorem 7. Let CI be an open Riemann surface with finite genus and at 
most countably many ends. Let Ki, i € / be some of the ends of ft, none 
of them punctures, where the index set L may be either finite or countably 
infinite. Suppose that for some e > and for each i £ J there is a ring 
domain Ai contained in Vt, with modulus at least e, that separates the end 
Ki from the other ends of CI. Then Q is conformally equivalent to a domain 
(p(fl) in a compact Riemann surface S, with complementary components Li 
corresponding to Ki, such that ip(£l)U\J ieI Li is a circle domain in S and for 
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each i £ I, L{ is either a singleton or a closed geometric disc with respect 
to the domain ip(tl) U Lj. Such a pair (S,tp) is unique up to conformal 
isomorphism. 

Proof. We begin by replacing CI by a conformally equivalent circle domain 
CI* in a compact Riemann surface R, as in Theorem [Bj Apply the construc- 
tion of section [5] to each of the components Ki, i £ I, to obtain a Beltrami 
differential fJ>{z)% supported on (J-fQ. In the case where I is infinite, the 
separation condition guarantees that \\fJ.\\oo < l f as it did for Theorem 
Solving the Beltrami equation fi(z)f z = on R gives a new compact Rie- 
mann surface S' with a homeomorphism / : R — > R' , conformal away from 
the preimages of the Ki, such that f{Kj) is a closed hyperbolic disc with 
respect to the hyperbolic metric of f(U) U f(Ki). The domain / (U U |J Ki) 
is countably connected so by Theorem [B] there is a conformal homeomor- 
phism g : f (U L)\J Ki) — > CI*, where CI* is a circle domain in a compact 
Riemann surface S. The image CI = g(f(U)) is the required domain, and 
Li = g(f(Ki)). Indeed, g is conformal on U U f(Ki), so it is an isometry 
from the hyperbolic metric of U U f(Kj) to the hyperbolic metric of CI U L,;. 

To prove uniqueness, suppose that / : CI — > CI' is a conformal homeo- 
morphism between two domains 0, C S and 0' C S' , each satisfying the 
conditions of the theorem. Let the complementary components of CI' be 
i £ I. Apply the argument of fj3]to extend / to a conformal homeomorphism 
of countably connected circle domains / : Q U [J Li — t CI' U [J L\. This must 
be the restriction of a conformal map from S to S' by the final part of He 
and Schramm's Theorem iBl □ 

7. Numerical approximation using circle packings 

The rest of this paper concerns the practical numerical approximation of 
finitely-connected intrinsic circle domains in the sphere. The idea is roughly 
to follow the steps of the existence proof, performing each step numerically. 

Firstly, we need a way to approximate circle domains having Fuchsian 
symmetry. This already presents us with a choice of methods, as there are 
several numerical methods for computing circle domains. Finitely connected 
circle domains can be approximated using iterated Riemann mapping, the 
individual steps of which can be carried out using a number of different 
methods for computational conformal mapping. However, we must deal 
with an infinitely connected circle domain, so it seems at first sight that we 
must make decisions about which complementary components to ignore at 
any particular stage of the approximation. We will see that this truncation 
can in fact be avoided. 

Secondly, if we follow the existence proof directly, it seems that we need 
to compute a quasiconformal extension numerically, find its dilatation and 
then solve a Beltrami equation on a larger domain. In practice, these steps 
can be combined into a single conformal welding step. 

A final practical difficulty may arise in converting the output of each step 
into the input for the next step, especially when the steps require different 
types of grid. 

Fortunately, circle packing is a numerical approximation method that of- 
fers solutions to all of these challenges. Although circle packing is not often 
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used for high-precision conformal mapping problems on account of its rela- 
tively slow convergence, it is particularly well suited to the approximation 
of circle domains, especially those with Fuchsian symmetry. It is also a 
practical tool for numerical conformal welding. We were able to carry out 
all of the steps of an approximation procedure for computing intrinsic circle 
domains using the CirclePack software [9] written by Ken Stephenson et 
al. For an introduction to circle packing and its use as a computational tool, 
we refer the reader to Stephenson's monograph [8]. 

7.1. Overview of circle packing. Let T be a finite graph embedded in 
the plane. A circle packing of T in the plane is a collection V of circular discs 
C v with disjoint interiors, one corresponding to each vertex v of T, such that 
whenever vertices v and w are adjacent in T the discs C v and C w are tangent. 
We call T the nerve of V . We may also circle pack in the sphere, using discs 
in the spherical metric, or in the hyperbolic plane, using hyperbolic circles 
and possibly also horocircles on the boundary. The carrier of V is the 
geometric complex formed by connecting the centers of neighboring circles 
by geodesic segments. According to the Koebe-Andreev-Thurston theorem, 
when T is any triangulation of the sphere, there exists a circle packing of 
T, and it is unique up to Mobius maps and reflection. In this case, the 
carrier of V is the entire sphere. Since we want to use circle packings to 
approximate conformal structures, we remove the reflection ambiguity by 
imposing a fixed orientation. 

There is a beautiful algorithm due to Bill Thurston for computing a circle 
packing from a given nerve. It works by removing one vertex and packing the 
remaining triangulation into the unit disc. This is achieved by computing 
the hyperbolic packing label for the packing; this is the function which assigns 
a hyperbolic radius to each vertex of the triangulation, in such a way that 
the boundary circles are given infinite radius, so that they correspond to 
horocircles, internally tangent to the unit circle, and all interior circles have 
finite hyperbolic radius. There is a unique packing label that results in an 
angle sum of 2n at each interior vertex. Thurston's algorithm approximates 
the correct label as the limit of a pointwise increasing sequence of labels. 
Once the correct packing label is known, the circles of the given radii can 
be laid out in the hyperbolic plane iteratively so that each satisfies the 
appropriate tangency conditions. 

7.2. Outline of the approximation algorithm. Here we outline a nu- 
merical method for approximating intrinsic circle domains. Let O be a given 
finitely connected domain in the Riemann sphere. Begin by taking a circle 
packing V such that the interiors of all the circles of V are contained in 
f2 and the complementary components of O are separated by the carrier of 
V . All subsequent calculations depend only on the nerve of V, so it is in 
this step that we have captured a discrete approximation to the conformal 
equivalence class of f2. The quality of the final approximation of the intrinsic 
circle domain £1* conformally equivalent to f2 will depend on the mesh of 
V (the size of its largest circle) and the maximum distance of dVt from the 
carrier of V . 



INTRINSIC CIRCLE DOMAINS 



15 



We perform a sequence of circle packing computations in order to con- 
struct a triangulation of the sphere together with an embedding of the nerve 
of V as a subcomplex. When we compute the circle packing Q of the spher- 
ical triangulation, the carrier of the embedded sub-packing gives an approx- 
imation to the intrinsic circle domain f2 conformally equivalent to U. We 
can interpolate the mapping of circle centers of V to the corresponding cir- 
cle centers in Q to give a polyhedral embedding of the carrier of V into the 
carrier of Q. Since we are mapping from the plane to the sphere, we cannot 
map in a piecewise affine fashion, so instead we map in a piecewise affine 
fashion to the polyhedron in R 3 whose vertices are the vertices of Q in S 2 , 
and then project radially outwards to S 2 . This map gives a homeomor- 
phism which approximates the conformal map from U to ft and is locally 
quasiconfor mal . 

We will explain informally why each step of the computation provides an 
arbitrarily good approximation to a conformal map occurring at the cor- 
responding step in the constructive proof of Theorem [TJ Note that the 
computation consists of finitely many circle packing steps. In each step we 
are approximating the solution of a problem whose solution depends contin- 
uously on its data. This should suffice to prove local uniform convergence 
of the approximations to the conformal map from £1 to £1*, but we will not 
attempt to estimate the rate of convergence. 

7.3. Approximating circle domains via circle packing. Suppose we 
are given a bounded finitely connected domain Q in the complex plane. We 
find a sequence of circle packings V n in f2 whose carriers exhaust £1. To 
construct V n we cut out a portion of the regular hexagonal circle packing 
with circles all of radius e = 2~ n . This means that we keep only those circles 
whose interiors are entirely contained in Q. Stephenson [8] describes this as 
using 80, as a 'cookie cutter'. We have to take some care at the boundary: we 
retain only the largest connected component of the hexagonal packing that 
remains, discarding any peripheral islands and iteratively removing vertices 
of degree one until we are left with a connected triangulation. When we 
take n large enough, the boundary components of this triangulation will be 
in one-to-one correspondence with the complementary components of Q. 

Next, to each cycle of boundary vertices of V n we adjoin a new 'ideal' 
vertex, adjacent to all of the vertices in the cycle, and in this way we obtain 
a triangulation of the sphere. By the Koebe-Andreev-Thurston theorem 
there is a circle packing whose nerve is this triangulation, and it is unique 
up to Mobius transformations. Consider what happens to this spherical 
packing as we let n tend to oo. Using quasiconformal distortion estimates 
for packings of bounded vertex degree, one can show that after a suitable 
Mobius normalization, the circles corresponding to the added ideal vertices 
converge to limit circles of positive radius, while the maximum radius of the 
all the other circles tends to zero uniformly. 

Let f n be the polyhedral map from the hexagonal cut-out packing V n to 
the carrier of the spherical packing, normalized so that three chosen points 
in are fixed by f n . It turns out that as n — > oo, f n converges locally 
uniformly on to a conformal map from 0, to a circle domain. To prove 
this, we can use the hex packing lemma of Rodin and Sullivan [6]. This says 
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that the quasiconformal dilatation of f n converges locally uniformly to on 
Cl. The Rodin-Sullivan bound for the quasiconformal dilation on a given 
face of the triangulation depends only on the number of layers of hexagonal 
packing surrounding it, and tends to as the number of layers tends to 
infinity. Let Cl' be any subdomain of Cl bounded by Jordan curves in Cl, 
such that Cl' is homeomorphic to Cl. The dilatation bound shows that the 
images of the maps /„ restricted to Cl' remain within a compact subset of the 
Teichmiiller space of multiply connected domains homeomorphic to Cl. Any 
subsequential limit of the sequence f n must be a conformal homeomorphism, 
and its image must be a circle domain. Given the normalization, there is 
a unique such homeomorphism. It follows that the sequence f n converges 
locally uniformly as n — > oo. 

We can augment the image packings by inverting across all the 'ideal' 
circles, adding circles of degree 4 in the resulting four-sided interstices to 
maintain a triangulation. By considering moduli of ring domains in the 
resulting packing it is possible to show that the maximum radius of any 
non- ideal circle in the image packing also converges to as n — > oo. 

We have now seen how to compute an arbitrarily good approximation to 
the circle domain Cl* conformally equivalent to Cl. We may assume that Cl* 
is contained in the plane. 



Figure 3. A circle packing of a circle domain in the Rie- 
mann sphere. Both the domain and the circle packing have 
sixfold symmetry. 



7.4. Circle domains with Fuchsian symmetry. Suppose B is one of 
the complementary components of Cl*. Let ir : D — > Cl* U B be an analytic 
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universal covering map. To compute an approximation to the circle domain 
conformally equivalent to H>\tt~ 1 (B), we use the same trick again. We have 
a triangulation of $7*, to which we add just one ideal vertex vb corresponding 
to the complementary component B. Now we apply Thurston's algorithm 
to compute the hyperbolic packing label for the resulting multiply connected 
triangulation. If we lay the circles out according to this packing label, we 
get a Fuchsian monodromy group. Call the resulting packing Vb- The 
circle corresponding to and its translates under the monodromy group 
approximate the bounded components of the required circle domain. For the 
conformal welding step that follows we will only need to know the hyperbolic 
radii for the circles corresponding to the vertex vb and its neighbors. It is 
not necessary to perform the layout routine. 




Figure 4. A circle packing approximation of the intrinsic 
circle domain conformally equivalent to the circle domain in 
Fig. [3l with the circle packing of the same complex shown. 

7.5. Conformal welding via circle packing. Given a Jordan curve 7 in 
the Riemann sphere, we can compute the Riemann mappings of the unit 
disc onto the interior and the exterior of 7. The boundary values of these 
Riemann mappings give two different continuous maps of the unit circle 
onto 7, (pi n t and 4> ex t- The shape of 7 is encapsulated in the boundary 
correspondence (j)~^ t o(j) int ^ which is an orientation-reversing homeomorphism 
of the unit circle. Conformal welding is the process of recovering 7 (up to a 
Mobius map) from the boundary correspondence. 

The CirclePack software includes a package for performing approximate 
conformal welding, following the procedure described by Williams [10J. The 
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idea is to paste two combinatorial closed discs together along their bound- 
aries, according to a best possible combinatorial approximation of the given 
boundary correspondence, to produce a triangulation of the sphere. We 
then compute the circle packing for this triangulation. For example, we 
may start with two maximal packings of the unit disc (meaning that the 
boundary circles are internally tangent to the unit circle), and use the arc 
length around the unit circle as a guide for pasting together the two trian- 
gulations, introducing new vertices where necessary to keep control of the 
vertex degrees. 

In our application to intrinsic circle domains, we weld a combinatorial 
closed disc to the boundary cycle corresponding to each complementary 
component B of f2, using the arc length around the circle corresponding to 
vb in the packing Vb and the trivial boundary correspondence e %d i— >■ e~ ld . 
A suitable choice of combinatorial closed disc would be a large section of the 
regular hexagonal packing, for then the hex packing lemma of Rodin and 
Sullivan can be applied to show that the desired welding and the computed 
welding differ by a quasiconformal map which has small dilatation except 
on an annulus of small modulus and small area covering the cycle of edges 
along which the welding is performed. After performing a welding for each 
boundary component, we simply apply Thurston's algorithm to compute a 
circle packing of the resulting spherical triangulation. In the case where 
the ambient Riemann surface R has positive genus, we would compute a 
packing of the universal cover of the welded complex at this stage to get a 
circle packing approximation to the Riemann surface S. 

The finished spherical circle packing gives an approximation to the map ip 
and to the desired conformal structure on the interior of each complementary 
component of <p(£l), and on the interior of the image of ip{Q). 

7.6. Convergence. We sketch a proof that the approximation scheme de- 
scribed above does converge locally uniformly on the domain Q* to the 
conformal mapping to an intrinsic circle domain, after suitable normaliza- 
tion. The strategy is to show that every map appearing in Figure [1] is well 
approximated by the corresponding polyhedral map between circle packings. 
Then since each step depends continuously on its input data, we will find 
that the polyhedral map from V* to S 2 converges locally uniformly on CI* 
to the desired conformal mapping. 

On the subcomplex of the packing representing the intrinsic circle domain, 
the hex-packing distortion estimate of Rodin and Sullivan shows that the 
polyhedral mapping is quasiconformal with dilatation that converges locally 
uniformly to zero. We can apply a similar estimate on the interior of each 
complementary component. 

It remains to show that the conformal welding step is well approximated 
by the discrete conformal welding. To do this we have to control the con- 
formal modulus of the image of a narrow annular neighborhood of each 
boundary component, where the circles are not deep enough in the hex 
packing for the Rodin-Sullivan estimate to give the bound we need. To do 
this it suffices to bound the maximum degree of the vertices appearing in 
the final triangulation. This is because the Rodin- Sullivan ring lemma allow 
the quasiconformal distortion of polyhedral maps between circle packings 
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to be bounded in terms of the maximum degree. We can arrange that all 
vertices have degree 6 except those involved in the combinatorial welding. 
The largest degree arising when we weld using the hyperbolic radii of the 
neighbors of vb in the packing Vb can be bounded in terms of the ratio 
between the radii of the largest and smallest circles adjacent to vb- Since 
we began with a hex-packing of with circles all of equal radius, this ra- 
tio is controlled in the limit by the maximum and minimum value of the 
derivative on the boundary dB of the conformal map h o ir" 1 . 




Figure 5. A circle packing of a circle domain with four com- 
plementary components. 

7.7. Practical considerations. Examining the figures, the reader will note 
that in our numerical experiments we have not restricted ourselves to using 
hex packings as our starting point. Instead we have used packings with few 
layers of degree 7 vertices near the boundary, so that there can are many 
very small circles on the boundary without needing a very large number of 
circles in total. This was done to reduce the errors from the combinatorial 
welding step. We pay for this by having larger circles far from the bound- 
ary, and losing the rigorous distortion bounds, but this seems a worthwhile 
trade since we expect the Schwarzian derivative of the map that we are 
approximating to be small on this region. In Figures [3] and H] there is a 
circle of degree 15 fixed by a rotational symmetry of the packing of order 3. 
The Schwarzian derivative of the map we are approximating must certainly 
vanish at the center of this circle, by symmetry considerations. Apart from 
two circles of this type, the maximum degree in the packing is 10. In the 
packings in Figures [5] and [6l the maximum degree is 7. 
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Figure 6. The spherical circle packing used to compute the 
intrinsic circle domain conformally equivalent to the circle 
domain shown in Fig. [5j 

Finally we comment that it may be more satisfactory from the point of 
view of computation, and also to simplify the convergence proof, to perform 
the discrete conformal welding using packings with specified overlap angles 
rather than restricting ourselves to tangency packings. However, we have 
not implemented this in practice. 
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